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SOLUTION OF THE POTENTIAL IN A SEMICONDUCTOR 
WITH EXPONENTIALLY DEPTH-DEPENDENT CONDUCTIVITY AND 
APPLICATION TO FOUR -POINT -PROBE MEASUREMENTS 

By R. K. Franks* and J. B. Robertson 
Langley Research Center 

SUMMARY 

An exact solution has been obtained for the potential due to a point current source 
on a semiconductor whose conductivity varies exponentially with depth. The solution was 
obtained by solving the continuity equation in cylindrical coordinates and applying the 
appropriate boundary conditions. This solution was then applied to the interpretation of 
the "four -point-probe” method of measuring semiconductor conductivity. 

INTRODUCTION 

The four -point probe is a favorite tool for measuring conductivity because it elim- 
inates the problem of contact resistance by use of null -measurement techniques. How- 
ever, previous theory has permitted the interpretation of measurements of bulk electrical 
conductivity by the four -point -probe method only for media of homogeneous conductivity 
(ref. 1). Some approximations have been made for nonhomogeneous conductivity such as 
small -amplitude sinusoidal variations (ref. 2). Extension of the theory to permit four- 
point -probe measurements of conductivity of nonhomogeneous media is very desirable, 
especially for persons studying diffusion of electrically active impurities. 

This paper presents a solution for the problem of a semiconductor of finite thick- 
ness whose conductivity depends exponentially upon depth. This solution may also serve 
as a basis for approximation in many of the actual situations involving depth -dependent 
conductivity. 

Included in the present paper is an appendix by Carl L. Fales, Jr., of the Langley 
Research Center, which explains the behavior of the electric field at small values of 
current -probe radius. 


*NASA-ASEE Summer Faculty Research Fellow from Applachian State University. 


SYMBOLS 


a exponential coefficient 

Ai,A 2 ,A 3 

b > b 1> b 2> b 3 

> constants 

C>Ci,C2,C3, ( - ; 4 
d 1> d 2> d 3> d 4 


E(r,z) 

E r (r,z) 

E z (r,z) 

F(a,S,T) 

h 

I 

*o 

Jo 

J(r,z) 

J r ( r ,z) 

Ko 

r 


base of natural system of logarithms, 2.718 

electric field vector 

r -component of electric field 

z -component of electric field 

correction factor 

depth of current probe 

current 

modified Bessel function of the first kind, zero order 
Bessel function of the first kind, zero order 
current density 

r -component of current density 

modified Bessel function of the second kind, zero order 
cylindrical coordinate 
radius of current probe 
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R solution to radius -dependent portion of equation for the potential 

5 probe spacing 

t time 

T thickness of semiconductor 

V potential difference between inner probes 

Y 0 Weber's Bessel function of the second kind, zero order 

z cylindrical coordinate 

z' point on vertical axis 

Z solution to depth -dependent portion of equation for the potential 

a 2 constant of separation 

6 Dirac delta function 

9 cylindrical coordinate 

p(r,z) free -charge density 

cr(z) conductivity 

<3>(r,z) coulombic potential 

if/ n set of eigenfunctions (n = 1, 2, . . °°) 

POTENTIAL OF A POINT CURRENT SOURCE 

The essential part of the problem is to determine the exact solution for the potential 
due to the presence of a point current source on the surface of a semiconductor whose 
conductivity depends exponentially upon depth; that is, 

ct(z) = cr(0)e az (1) 
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Cylindrical coordinates (r, 9, z) are used and the semiconductor is bound at the planes 
z = 0 and z = T and at the cylinder r = r c . (See fig. 1.) 


r = 0 



Figure 1.- Semiconductor with current probe inserted distance z = h into 

cylindrical cavity. 

The cylindrical current probe is inserted a distance h into the cavity as shown in 
figure 1. With the assumption of uniform current density, the boundary conditions on the 
cylindrical surface are 


E r (r c ,z) = 


I 

27rr c ho-(z) 


(0 < z < h) 


E r( r o z ) = 0 


(z ^ h) 


( 2 ) 


where E r (r c ,z) is the radial component of the electric field at the surface r c , I is 
the current carried by the probe, and cr(z) is the conductivity. 

The radial component of the electric field must go to zero as r approaches infinity; 
that is, 


E r (°°,z) = 0 


(3) 
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The z -component of the electric field must be zero on the planes z = 0 and z = T; 
that is, 

E z (r,0) = E z (r,T) = 0 (4) 

In the limit, as both r c and h approach zero, a situation arises which is equiv- 
alent to a point current source on the surface of the semiconductor. It has been assumed 
that if minority carriers are injected, they will recombine close enough to the electrode 
so that the conductivity will not be affected. 

The continuity equation now applied is 

V • J(r,z) = V • |a(z)E(r,z)] = - -^-p(r,z) (5) 

where J(r,z) is the current density and p(r,z) is the free -charge density. By using 
equation (1) for the conductivity and expressing the electric field as the negative gradient 
of the potential (E = -V$(r,z)), equation (5) under time -independent conditions becomes 


V 2 $(r,z) + a 4>(r,z) = 0 


( 6 ) 


Assumption of a product solution of equation (6) of the form <i>(r,z) = R(r)Z(z) and 
separation of variables yield 


r 2 -d__ R + r iL R - a>2 r 2R = o 
dr 2 dr 


(7) 


and 


-d— Z + a -d- Z + a 2 Z = 0 


( 8 ) 


dz 2 dz 

where a 2 is the constant of separation. 

There are three physically possible solutions to equation (7), each corresponding to 
a different value of a 2 . These solutions (ref. 3, sec. 4.8) are 


R 1 = A iI 0 ( Q!r ) + B l K 0 (ar) 

R 2 = A 2 + B 2 r 

r 3 = A 3 J 0 (o!r) + B 3 Y 0 (o!r) 


(at 2 > 0) 
(a 2 = 0) 
(of 2 < 0) 


"\ 


(9) 


J 


where I Q , Kq, J q , and Y 0 are zero -order Bessel functions. The boundary condi- 
tion E r (°°,z) = 0 requires that Aj = 0. 
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There are four possible solutions to equation (8) corresponding to different values 
of q» 2. These solutions are 


Zj = e -az /^(Ci sin kz +Dj cos kz) 


a 2 > ^ 
4 


(10a) 


where k = |/q ! 2 - 


Z 2 = C 2 e" az / 2 + D 2 ze _aZ//2 

II 

CM 

Z 3 = e- az / 2 (c 3 eg z +D 3 e-g z ) 

^^>a 2 >oj; (a 2 < 0) 

where g = - cn 2 , and 



(10b) 

(10c) 


Z 4 = C 4 + D 4 e" az 


(a 2 = 0) (lOd) 


The boundary condition E z (r,0) = E z (r,T) = 0 requires that C2, C3, D2, D3, and D 4 

be zero, that C 1 = ^ D l, and that k = M where n = 1, 2, . . .,00. Therefore, 


Zi,„ - e ' aZ/2D n(i£ sin f z + cos f z) 

z 2 = 0 
z 3 = 0 
z 4 = c 4 




( 11 ) 


The solution of equation (6) is a linear combination of all product solutions 

OO 

*(r,z) = £ A n e- az/2 K 0 (Q! n r)^ sir 


n=l 


sin z + cos z) + B In r + C (12) 


where 


“n = faf) + ^ = ^jfnV+aV 
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By using boundary conditions at r = r c and the orthogonality of the terms in 
E r (r c ,z) with respect to the weighting function e az over the interval z = 0 to z = T, 
it can be shown (ref. 3, sec. 5.6) that in the limit of small r c and h 


and 


B = 


al 

2mx(0)(l - e a ^) 


(13) 


A„ = 


4 n 2 7r 2i 


7TCT (0)T(a 2 T 2 + 4n 2 5T 2 ) 

Substituting for A n and B in equation (12) gives 


(n = 1, 2, . . oo) 


(14) 


4>(r,z) = 


Ie 


-az/2 


7ra(0)T 


I ^7^2 •*» ¥ 


nir 

z + cos z 
T 


al In r 


2ttct( 0)(1- e aT ) 


+ C 


(15) 


Equation (15) is the solution which has been sought. It is an equation for the poten- 
tial in cylindrical coordinates, independent of 9 and in terms of the probe current, 
semiconductor thickness, and variation in conductivity. The behavior of the electric 
field, as determined from this expression for the potential, in the limit of small r, is 
treated in the appendix. 

This solution for the potential inside a semiconductor should be of interest to 
designers of semiconductor devices since it describes the field variations within the bulk 
material. The authors' prime interest in this solution is its use in interpreting four- 
point -probe measurements. 

APPLICATION TO FOUR-POINT PROBE 

The present solution (eq. (15)) is now applied to the case of a four -point probe on a 
flat semiconductor. The common four -point probe consists of four colinear, equispaced 
point probes (fig. 2), The present treatment is justified if the probe contact radius is 
small compared with the probe spacing and if the lateral boundary distance is large com- 
pared with the probe spacing. 
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Fig-ure 2.- Schematic diagram of four-point probe and associated circuit in 
contact with semiconductor surface. 

If probe 1 introduces current +1 and if probe 4 introduces current -I, then the 
potential difference between probes 2 and 3 due to the current source at probe 1 will be 

V + = 4>(S,0) - $(2S,0) 

and the potential difference between probes 2 and 3 due to the current source at probe 4 
will be 

V = -|$(2S,0) - $(8,0)] 

The total potential difference V between the two inner probes will then be 

V = V + + V" = 2 ji(S,0) - 4>(2S,0)| (16) 
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Use of the present solution for the potential (eq. (15)) gives 


r « 


V - -2L./1 Y 4n 2 ?r 2 K /_S_l/ 4n 2„2 . a 2 T 2 
' ™«» ^ »4 4n V + a*T* r°\2T ^ + * T 


K 0 (S.^4n 2 jr 2 + a 2 T 2 


a In i 


/J 2(l-e aT )J 


( 17 ) 


which can be solved for the conductivity at the surface cr(0). The resulting equation is 

a(0) = zkv F ( a > s > T ) 


where 


F(a,S,T) 



4n 2 7r 2 

4n 2 7r 2 + a 2 T 2 



+ a 2 T 2 


- K 0 ^4n 2 ?r 2 + a 2 T 2 j 


2aS In | 

+ rr^ 


(18) 


The factor F(a,S,T) is a correction factor taking into account the semiconductor thick- 
ness, the probe spacing, and the variation of the conductivity. 

The equation for conductivity (eq. (18)) reduces, appropriately, to the equation for 
conductivity in a sheet conductor (ref. 4) for very small semiconductor thickness T 
and reduces to the equation for homogeneous bulk conductivity (ref. 1) when the coeffi- 
cient a goes to zero. The summation term in F(a,S,T) converges rapidly. 

Figure 3 presents F(a,S,T) as a function of semiconductor thickness T in terms 
of probe spacing S and for various values of the exponential coefficient a. 
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Figure 3-- Correction factor F(a,S,T) as a function of semiconductor thickness T in 
terms of probe spacing S and for various values of exponential coefficient a. 


CONCLUDING REMARKS 


The solution for the potential in a semiconductor with exponentially depth dependent 
conductivity allows for the first exact interpretation of four -point -probe resistivity mea- 
surements on semiconductors of nonhomogeneous conductivity. Although a semiconductor 
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with a true exponential conductivity dependence may not be encountered, the treatment 
presented herein serves as a basis for approximation for other depth dependent conduc- 
tivities such as those resulting from diffusion of impurities. This solution for the poten- 
tial inside a semiconductor should also be of interest to designers of semiconductor 
devices since it describes the field variations within the bulk material. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. September 1, 1971. 
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APPENDIX 


ELECTRIC FIELD AT SMALL VALUES OF r c 

By Carl L. Fales, Jr. 

Langley Research Center 


A comment about the behavior of the r -component of the electric field in the limit 
of small r seems desirable. 

In the limits h — 0 and r c — 0, 




_ I 


27tct(0) 




n=l 


n-faT') 

V2nW 


2l _1 


-az/2[~ d_ 
" dr 


Ko(«r)j^ 


aT 

2n7r 


sin ~~ z 


+ cos 


¥■) 



For r finite, there exists a number N such that Q! n r » 1 for all n > N. 
Therefore, 


KoJar) « Constant x ■■ e a?nI ' 
dr |/c^r 


and the series surely converges. 

On the other hand, it is expected that E r (0,z) = 0 for z > 0 by symmetry. How- 
ever, in the limit of small r, K^ar) — — and hence 


dr 


E, 


2ncr(0) 


r oo 


_ 

ly 2 

1 + 

f — ) 2 

1 A T 

\2n7r/ _ 

ln=l 




-az/2 / jT_ s ^ n M z + cos ^ z) + — — 
\2n7r T T ) e aT _ l 


= l I 

2ttct( 0) r 


X KM 


n=l 


Clearly since b n M -b 0 as n — °°, this series does not converge, at least in the sense 
of numerically summing a finite number of terms. 


12 


APPENDIX - Continued 


j2 j 

From equation (8), the operator L = -=-5 + a ^ can be shown to be Hermitian (due 


dz £ 


to the self -adjointness of L and the homogeneous boundary conditions on E r (r,z)). As 
a result, the set of eigenfunctions i// n (z) of the eigenvalue equation 


Li// n (z) = -a 2 ^ n 


(n = 1, 2, 


°) 


are orthonormal with respect to the weighting function e az and form a complete set 
(for functions obeying the same homogeneous boundary conditions). Hence, the repre- 
sentation for the Dirac delta function can be written as 


6(z - z') = Yj eaZ ' ’M 2 ) ’M 2 ’) 
n=l 


(*n real ) 


or 


6(z) = Y >M 2 ) 

n=l 


The properly normalized eigenfunctions are 





1/2 


and 


1/2 




1 * m 


e -az/2^^X. s i n M. z + cos M. 


Therefore, 


wj 

6 (z) - £ | 


n=l 




\2n7r/ 


T 1 

’ e ' aZ/2 (S sin f z + cos z) + 


and 


lim E r (r,z) = 
r— 0 


I 

2ircr(0) 


6(z) 
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APPENDIX - Concluded 


Thus, E r (0,z) = 0 for z > 0 as had been expected on physical grounds. Also, the 
component of current density 

Jr(r ’ z) = 2iF 6(z) 

gives the correct total current 


r™ J 27rr dz J r(r,z) = I 
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